Stability and chaos in real polynomial maps
Thesis presented by
Fermín Franco-Medrano, B.Sc.
As a requirement to obtain the degree of
Master of Science with specialty in Applied Mathematics
conferred by
Centro de Investigación en Matemáticas, A.C.
Thesis Advisor:
Francisco Javier Solís Lozano, Ph.D.
Guanajuato, GT, United Mexican States, June 2013.
Stability and chaos in real polynomial maps
Thesis presented by
Fermín Franco-Medrano, B.Sc.
As a requirement to obtain the degree of
Master of Science with specialty in Applied Mathematics
conferred by
Centro de Investigación en Matemáticas, A.C.
Approved by the Thesis Advisor:
Francisco Javier Solís Lozano, Ph.D.
Guanajuato, GT, United Mexican States, June 2013.
Abstract
We extend and improve the existing characterization of the dynamics of general quadratic real polynomial maps with coefficients that depend on a single parameter , and generalize this characterization to cubic real polynomial maps, in a consistent theory that is further generalized to real -th degree real polynomial maps. In essence, we give conditions for the stability of the fixed points of any real polynomial map with real fixed points. In order to do this, we have introduced the concept of Canonical Polynomial Maps which are topologically conjugate to any polynomial map of the same degree with real fixed points. The stability of the fixed points of canonical polynomial maps has been found to depend solely on a special function termed Product Distance Function for a given fixed point. The values of this product distance determine the stability of the fixed point in question, when it bifurcates, and even when chaos arises, as it passes through what we have termed stability bands. The exact boundary values of these stability bands are yet to be calculated for regions of type greater than one for polynomials of degree higher than three.
Acknowledgements
I have both the duty and pleasure of thanking the support of many people and institutions that helped me in the process of conducting my graduate studies. In the first place, I greatly thank the support of the National Council of Science and Technology of the United Mexican States (CONACYT) for their financial support as a national scholar during my studies. I thank my teachers and the personnel of the Center for Mathematical Research (CIMAT) in Guanajuato, for their patient and committed work from which I nurtured and which allowed me to grow in this Mexican institution. In particular, I thank the orientation and guidance of Professor Francisco Javier Solís Lozano, my thesis advisor; as well as Professors Lázaro Raúl Felipe Parada, Mónica Moreno Rocha and Daniel Olmos Liceaga, who contributed greatly to the improvement of this work with their comments and advice; and Dr. Marcos Aurelio Capistrán Ocampo for his support and orientation. I also thank Mr. José Guadalupe Castro López and my friends and classmates for their support and invaluable company. To my family and particularly to my mother, I must extend the greatest gratitude of all.
Chapter 1 Introduction
In this chapter, we will briefly address the basic concepts of the theory of discrete dynamical systems, so as the concept of chaos that can occur in these systems, as much as necessary to define the relevant problem and the goals of this thesis. Important references in this subject are [Devaney(1989), Holmgren(1994), Elaydi(2000)], among others, and all of the proofs omitted here can be found therein. This chapter is included with the aim of making the work self-contained, it includes only what is now the standard basic theory of discrete dynamical systems, and may be skipped if the reader is familiar with it.
This chapter has the following objectives: (i) familiarize with the fundamentals of the theory of discrete dynamical systems, particularly chaotic ones; (ii) show the importance of the parametric dependency in families of functions, particularly bifurcation theory; and (iii) present the Feigenbaum sequence and its universality property for unimodal maps as a basis to understand how chaos may arise through period doubling bifurcations.
1.1 Elementary concepts
The main objective of the theory of discrete dynamical systems is to understand the final or asymptotic behavior of an iterative process. If the process is continuous, a differential equation with time being the independent variable, then the theory intends to predict the ultimate behavior of the solutions in the distant future () or the distant past (). If the process is discrete, such as the iteration of a real (or complex) function , then the theory hopes to be able to understand the eventual behavior of the points when is large, where denotes the -th iteration (i.e., composition of the function with itself). That is, the theory of dynamical systems asks the question: where do points go and what do they do when they get there? In this work, we aim to answer a part of this question for one of the simplest kinds of dynamical systems: functions of a single real variable; in particular, for real polynomial functions, having as a precedent the work of [Solís and Jódar(2004)] for quadratic maps.
In the following, let , and .
Definition 1.1 (Orbits).
The forward orbit of is the set of points and is denoted by . If is a homeomorphism, we can define the complete orbit of , , as the set of points for , and the backward orbit of , , as the set of points .
Therefore, we can restate the basic goal of discrete dynamical systems as to know and understand all orbits of a map. Orbits can be quite complicated sets, even for simple non-linear maps. However, there are some specially simple orbits which play a central role in the dynamics of the whole system.
Definition 1.2 (Periodic Points).
The point is a fixed point of if . The point is a periodic point of period if . The smallest positive for which is called the prime period of . We denote the set of periodic points of period (not necessarily prime) by , and the set of fixed points by . The set of all iterates of a periodic point forms a periodic orbit.
An important question is the “range of influence” of a fixed point. That is to ask, what set of points, if any, will approach or tend to a fixed (or periodic) point under the iteration of .
Definition 1.3 (Stable Set or Basin of Attraction).
Let be periodic of period . A point is forward asymptotic to if . The stable set (or basin of attraction) of , denoted by , consists of all points forward asymptotic to .
If is not periodic, we can still define its forward asymptotic points requiring that when . Moreover, if is invertible, we can consider backward asymptotic points by taking in the last definition. The set of backward asymptotic points to is called the unstable set of and is denoted by .
Another useful concept comes from elementary calculus.
Definition 1.4 (Critical Points).
A point is a critical point of if . The critical point is non degenerate if . The critical point is degenerate if .
Critical points will prove later to be important in relation to the stable set of a fixed point, and to determine how many attracting orbits a map may have.
Degenerate critical points can be either maxima, minima or saddle points. Non degenerate critical points can only be maxima or minima. A diffeomorphism cannot have critical points, but its existence in non-invertible maps is one of the main reasons why this type of maps are more complicated.
The goal of the theory of discrete dynamical systems is then to understand the nature of all orbits and to identify which are periodic, asymptotic, etc. In general, this is an impossible task. For example, if is a quadratic polynomial, then explicitly finding the periodic points of period is equivalent to solving the equation , which is a polynomial equation of degree . Computer calculations do not help in this case either, since rounding errors tend to accumulate and make many periodic points “invisible” to the computer. Therefore, qualitative and geometrical techniques must be employed to understand the nature of a system.
1.2 Hyperbolicity and criteria for stability
A key issue once we find periodic points, is to determine the behavior of orbits near these points; such investigation is what is called stability theory. Once this is done, we can state whether a periodic point is attracting or repelling (or neither). To determine this, we can part from the following
Definition 1.5 (Stability).
Let be a fixed point of . Then,
-
1.
is said to be stable if for any there exists a such that if then for all integers . Otherwise, the fixed point is called unstable.
-
2.
is called attracting if there is a such that if then .
-
3.
is asymptotically stable if it is both stable and attracting.
Therefore, intuitively, a fixed point is stable if nearby orbits remain near the fixed point, and you can always find an interval around the fixed point such that orbits starting within the interval will remain arbitrarily close to the fixed point. The opposite is that you cannot do that; i.e. there is an interval around the fixed point such that no matter how close to the fixed point you start the orbit from, the orbit will go outside the interval. The case of an asymptotically fixed point is particularly important, since in this case all the orbits near a fixed point will approach it in the limit. The next step is to be able to determine when a fixed point is stable or unstable. To do this, we will need the
Definition 1.6 (Hyperbolic Point).
Let be a periodic point of period . The point is hyperbolic if . Otherwise, it is called a nonhyperbolic fixed point. The number is called the multiplier of the periodic point.
Example 1.1.
Consider the diffeomorphism . There are three fixed points: and -1. Note that and . Therefore, each point is hyperbolic.
As it turns out, hyperbolic points are easy to understand in terms of whether they are attracting or repelling, as they are always one of the two.
Definition 1.7 (Attractor and repellor).
Let be a hyperbolic point of period . If , the point is asymptotically stable and is called an attracting periodic point or more shortly, an attractor; if , is unstable and it is called a repelling periodic point or, simply, a repellor.
Occasionally, the terms sink and source may also be used to refer to attractors and repellors, respectively, stemming from the continuous dynamical systems analog. Now, if a periodic point is attracting, it must have an interval around it in which it is so; this is stated in the following
Theorem 1.1.
Let be a hyperbolic fixed point of . Then if is an attractor, there exists an open set , with ; on the other hand, if is a repellor, there is an open set , with .
This theorem explains the choice of the terms “attracting” an “repelling”, since attractors and repellors have around them a basin of attraction and an unstable set, respectively.
It is straightforward to see that is an invariant set under the action of the map. A similar result is true for periodic points of period . In this case, we have an open interval around that is mapped into itself by .
1.3 The Schwarzian derivative
Here, we will discuss a mathematical tool that will turn out useful in determining some properties of maps in terms of its periodic points structure.
Definition 1.8 (Schwarzian Derivative).
The Schwarzian derivative of a function at the point is given by
In particular, it was first used by [Singer(1978)] in [Singer(1978)] to address the question of how many attracting periodic points a differentiable map can have [Singer(1978)] (see Singer’s Theorem below). It can also be used to determine the nature of non-hyperbolic periodic points [Elaydi(2000)]. Mappings with negative Schwarzian derivative present particular dynamical properties that will interest us in this work.
Theorem 1.2.
Let be a real polynomial. If all roots of are real and distinct, then for all .
The following theorem is useful when analyzing topologically conjugate maps (a concept to be defined further below in 1.9), which in turn is useful to determine the stability properties of one map in terms of the known such properties of another map (the conjugate map).
Theorem 1.3.
Suppose and . Then .
Also, this tells us that iteration of a map does not change the sign of its Schwarzian derivative, i.e. the property is preserved. This immediately takes us to
Corollary.
Let . Then for all .
The assumption of has several surprising implications for the dynamics of a one-dimensional map, as we will see below. Another useful result dealing with the geometry of a map is
Lemma 1.4.
If , then cannot have a positive local minimum or a negative local maximum.
The proof of this last lemma is straightforward from the definition of the Schwarzian derivative and elementary calculus. Note that this lemma tells us that, if , between any two consecutive critical points of , its graph must cross the -axis; this in turn means that there must be a critical point for between these two points, i.e. there must be a maximum or minimum of between any pair of its inflection points.
The Schwarzian derivative can also be used to determine whether a fixed point attracts a critical point of a map.
Lemma 1.5.
Let and be fixed points of a continuously differentiable map with and such that on the open interval . If , then has a critical point in .
The last lemma is used in the proof of the following important result due to [Singer(1978)] ([Singer(1978)]).
Theorem 1.6 (Singer’s Theorem).
Suppose ( is allowed). Suppose has critical points. Then has at most attracting periodic orbits.
This last theorem allows to find directly from an upper bound for the number of attracting periodic points. This differs from the information provided in the analysis by Sarkovskii’s theorem (see section 1.5 below) in that the latter tells us how many periodic points there are, regardless of whether they are attracting or not and, also, we have first to find some periodic point and its period to be able to apply Sarkovskii’s theorem.
1.4 Nonhyperbolic fixed points
The stability criteria for nonhyperbolic fixed points involve the Schwarzian derivative defined in 1.3. We will analyze the cases of the multiplier being equal to 1 and -1 separately. The unstated proofs of the following theorems can be found in the book by [Elaydi(2000)].
Theorem 1.7.
Let be a fixed point of a map such that . Then, if and is continuous, the following statements hold
-
1.
If , then is unstable.
-
2.
If and , then is unstable.
-
3.
If and , then is asymptotically stable.
Similarly, we have
Theorem 1.8.
Let be fixed point of a map such that . If is continuous, then the following statements hold:
-
1.
If , then is an asymptotically stable.
-
2.
If , then is unstable.
It is worth noting that a repellor is an unstable periodic point but not all such points are repellors since, a nonhyperbolic periodic point may also be unstable; e.g. a “semistable” nonhyperbolic fixed point is “attracting from the right (or left)” but “repelling from the left (or right, respectively)”.
When dealing with nonhyperbolic points, particularly when they are unstable, it is useful to account for the concept of “semistability”, separating stability “from the right” and “from the left”.
Definition 1.9 (Semistability).
A fixed point of a map is said to be semistable from the right (respectively, from the left) if for any there exists such that if (respectively, ) then for all . Moreover, if whenever (respectively, ) for some , then is said to be semiasymptotically stable from the right (respectively, from the left).
It is then easy to prove that
Theorem 1.9 (Asymptotical semistability).
Let be a map, one of its fixed points and with and . Then is
-
1.
Semiasymptotically stable from the right if ;
-
2.
Semiasymptotically stable from the left if .
1.5 Sarkovskii’s theorem
This theorem is of great beauty, as it constructs a new –apparently artificial– ordering of the natural numbers in order to determine what types of orbits may a map have. Consider then the following ordering of the natural numbers:
That is: first all odd numbers, then two times each odd number, times each odd number, times each odd number, etc., and finally, the powers of two in descending order. This is Sarkovskii’s ordering of the natural numbers and it easy to see it is exhaustive .
Having defined the above, we can state
Theorem 1.10 (Sarkovskii).
Let continuous. Suppose that has a periodic point of prime period . If in Sarkovskii’s ordering, then also has a periodic point of period .
Thus, Sarkovskii’s theorem gives a nearly complete answer to the question of how many periodic points and of which periods can a map have, although it does not tell us their nature in terms of their stability, nor does it tell us how to find such points.
Some important consequences of Sarkovskii’s Theorem can be summarized as follows:
-
•
If has periodic points with periods different from a power of 2, then has infinite periodic points of infinitely different periods.
-
•
If only has a finite number of periodic points, then necessarily all points have periods of powers of 2.
-
•
Period 3 is the “greatest” period in Sarkovskii’s ordering, therefore, it implies the existence of all other periods.
Sarkovskii’s theorem is without doubt a central result for the theory of discrete dynamical systems.
1.6 Bifurcation theory
We will now consider families of maps instead of single maps, as this is often useful in applications and is also the main object of study of this work. For example, in practice, the coefficients of a polynomial map that models some natural, physical, economic or social discrete process can only be determined up to a certain precision, which gives rise to a family of maps as a set bounded by the uncertainties of the coefficients.
The aim of bifurcation theory is then to study the changes that occur in a map when its parameters are varied. This changes commonly refer qualitatively to its periodic points structure and stability properties, but they can also involve other changes. Consider families of functions in one real variable which depend upon a single real parameter (uniparametric). Consider the bivariate function , where for fixed , . We will assume also that depends smoothly on . Some examples of such families of functions are:
-
•
.
-
•
.
-
•
.
-
•
.
With this in mind, we will next give a rather informal definition of bifurcation.
Definition 1.10 (Bifurcation).
Let be a uniparametric family of functions. Then, we say there is a bifurcation at if there exists such that if and satisfy and , then the “dynamics” of is different from the “dynamics” of . That is, the “dynamics” of the function changes when passes through the value .
The “changes in the dynamics” of the above definition refer to qualitative changes in the structure of periodic points of . To make the above notion precise, we must consider separately the different types of bifurcations that may occur. For real maps, the saddle-node and period doubling bifurcations are the most common. We will next discuss the concept of saddle-node bifurcation through an example and the period doubling bifurcation is discussed in more detail below.
Example 1.2 (Saddle-node or Tangent Bifurcation).
Let’s consider the family of functions , with . A bifurcation occurs when . We have three cases:
-
1.
If the graph of and do not intersect each other, so that has no fixed points. .
-
2.
When , the graph of tangentially touches the identity diagonal at . is the only fixed point. If , and if , .
-
3.
When , the graph of crosses the identity diagonal two times, at with and at with , with which we see they are an attractor and a repellor, respectively.
Figure 1.1 presents some of these characteristics graphically, and the corresponding phase portraits are presented in figure 1.2. Figure 1.3 presents a common and very useful way of depicting bifurcations of real maps, which will be frequently used in this work: a bifurcation diagram; in it, the horizontal axis represents the values of the parameter in which the family is varied, and the vertical axis represents the values of the periodic points present for each .
The above example suggests that bifurcations may occur when non-hyperbolic fixed points are present. Indeed, this the only case when bifurcations of fixed points occur (in the sense of one fixed point giving birth to two others), as the next theorem demonstrates.
Theorem 1.11.
Let be a one-parameter family of functions and suppose that and . Then there are intervals about and about and a smooth function such that and . Moreover, has no other fixed points in .
Proof.
Consider the function defined by . By hypothesis, and
| (1.1) |
using the other hypothesis. Then, by the Implicit Function Theorem [Marsden(1974)], there are intervals about and about , and a smooth function such that and for all ; that is, is fixed by for . Moreover, unless . ∎
The above theorem obviously holds for periodic points by replacing with .
1.7 Chaos
In this section we will precise what we understand by the mathematical concept of chaos. It must be said though, that there is not a single universally accepted definition for this concept, as variations occur; however, they preserve the same general “spirit”. We will here adopt one of the most popular definitions, due to [Devaney(1989)].
Before we can properly define the concept of chaos, we need some other auxiliary concepts.
Definition 1.11 (Topological transitivity).
Let be a metric space and a map. Then is said to be topologically transitive if for any pair of nonempty open sets there exists such that .
In plain words, what this definition states is that a map is topologically transitive if any arbitrarily small set of the domain of contains points that are mapped by to any other arbitrarily small set of the codomain in a finite number of iterations. An immediate consequence of topological transitivity is that the domain cannot be divided into two nonempty disjoint sets which are invariant under . It is easy to see that if has a dense orbit in then it is topologically transitive (the converse is also true for compact subsets of or ).
Another necessary concept for the definition of chaos is the following.
Definition 1.12 (Sensitive Dependence to Initial Conditions).
Let be a metric space and a map. is said to posses sensitive dependence to initial conditions if there exists such that for any and any neighborhood of , there exists and such that .
Intuitively, sensitive dependence to initial conditions means that no matter how close two points are to each other, they will eventually be separated by at least some distance under a finite number of iterations of .
With these preliminary definitions we can now state the main definition of this section.
Definition 1.13.
Let be a metric space and a map. Then is said to be chaotic on if
-
1.
is topologically transitive,
-
2.
the set of periodic points of is dense in , and
-
3.
has sensitive dependence on initial conditions.
So a chaotic map “mixes up” its domain set by means of its topological transitivity but everywhere you look there are periodic points arbitrarily close and, finally, it is practically impossible to predict orbits numerically; these are the three elements that make up a chaotic map, under this definition. Actually, it was later shown [Banks et al.(1992)Banks, Brooks, Cairns, Davis, and Stacey] that topological transitivity and denseness of periodic points imply sensitive dependence to initial conditions, but no other pair of conditions imply the other. However, for continuous maps on intervals of , topological transitivity implies that the set of periodic points is dense [Vellekoop and Berglund(1994)], so that in this case, by the above, topological transitivity alone actually implies chaos.
1.8 The period doubling route to chaos
Perhaps one of the most important types of bifurcations is the period doubling bifurcation. This type of bifurcation will play a central role in the main results of this work. The most commonly know map which presents period doubling bifurcation is the famous logistic map. The issue is that the logistic quadratic map has simple dynamics for but it becomes chaotic when (we will precise the notion of chaos in the next section). The natural question here is: how exactly is it that becomes chaotic as increases? Where do the infinitely many periodic points come from when is large? In this section we will see that “period doubling” bifurcation is central to this.
Sarkovskii’s theorem provides a partial answer to the question of how many periodic points emerge when the parameter is varied. Before can have infinite periodic points of distinct periods, it must first have periodic points of all the periods of the form with . Local bifurcation theory provides two typical ways in which these points may emerge: by means of saddle-node bifurcations or through period doublings. The question is then which type of bifurcations occur as becomes chaotic.
The typical way in which becomes chaotic is that undergoes a series of period doubling bifurcations. This is not always the case, but it is a typical “route to chaos”.
Theorem 1.12 (Period Doubling Bifurcation).
Suppose that
-
1.
for all in an interval around .
-
2.
.
-
3.
Then there exists an interval around and a function such that
but
Proof.
Let . Then, from hypotheses (1) and (2),
so that we cannot apply the implicit function theorem directly. Thus, to rectify the situation, we set
| (1.2) |
Straightforward calculation of the limits of the derivatives allow us to verify that is smooth and satisfies
| (1.3) | ||||
We now observe that, once again using hypotheses (1) and (2),
| (1.4) | ||||
Moreover, now using hypothesis (3),
| (1.5) | ||||
Therefore, by the implicit function theorem [Marsden(1974)], there exists a map defined on an interval around such that and . So that, for ,
Consequently, and then and thus is of period 2 for and for in an interval around . Notice that is not fixed by because we can apply theorem 1.11 in this case and, then, the curve of fixed points is unique.
On the other hand, from the chain rule we have
From which,
| (1.6) |
We know that, for , the denominator in (1.6) is different from zero and, on the other hand, the numerator is equal to and
| (1.7) | ||||
Where we have once again used hypothesis (2). Ergo, . ∎
Also, differentiating equation (1.6) again and evaluating in , it is easy to see that
| (1.8) |
We already know that the denominator in (1.8) above is not null and, also,
| (1.9) | ||||
where is the Schwarzian derivative of and, yet once again, we used hypothesis (2) from the theorem above. Therefore, substituting in equation (1.8) we have
| (1.10) |
So that, if then and the curve is either concave left or right. Along with this gives us an idea of the geometry of the period 2 curve near the bifurcation point . Notice that does not bifurcate in the sense of giving birth to another fixed point (as was the case in the pitchfork bifurcation) but, rather, continues to be the only fixed point in a interval around (see theorem 1.11).
In the case of the logistic map, the period doubling bifurcation repeats over and over again for the fixed point as the value of the parameter is increased, starting from :
-
•
Between , has a single stable fixed point, (of period ).
-
•
For , loses its stability and ’gives rise’ to an attracting cycle of period two (), whose orbital points separate from as increases.
-
•
For [Elaydi(2000)], the attracting cycle of period two loses in turn its stability and gives rise to an attracting cycle of period .
-
•
For , the period four cycle also loses its stability and bifurcates into an asymptotically stable cycle of period .
The above mentioned process repeats itself indefinitely and produces a sequence . The Feigenbaum sequence describes the values of the parameter for which period doubling bifurcations arise in a unimodal map [Feigenbaum(1978)]. The sequence converges exponentially. The important point to remark here is that the rate of convergence of these sequences (whose values depend on the specific unimodal map), is constant for every unimodal function and is given by , which gives rise to the so-called property of universality. In fact, its importance deserved naming it Feigenbaum’s constant, after [Feigenbaum(1978)] discovered it in [Feigenbaum(1978)] (though amazingly ignored for over a decade).
1.9 Topological conjugacy
Here we will define the important concept of topological conjugacy which relates the dynamics of two different maps by means of a third one. The importance of conjugacy is that it preserves the stability and chaotic properties of a map, so that one map can be analyzed by means of analyzing another one, presumably simpler, so that the analysis becomes easier.
Definition 1.14 (Topological Conjugacy).
Let and be two maps. Then and are said to be topologically conjugate, denoted , if there exists a homeomorphism such that . The homeomorphism is called a topological conjugacy. We also say that is -conjugate to to emphasize the importance of the homeomorphism .
It is easy to prove that conjugacy is an equivalence relation. Topologically conjugate maps have equivalent dynamics. For example, if and is a fixed point of then is a fixed point of ; and, also, for all . Similarly, periodic orbits of are mapped to analogue periodic orbits of .
Theorem 1.13.
Let be -conjugate to . If is chaotic on , then is chaotic on .
Theorem 1.13 will play a key role in the importance of the study of “canonical polynomial maps” in chapters 3, 4 and 5. The proof of this theorem can be found in [Elaydi(2000)], where it can be seen that, in fact, the topological conjugacy does not need to be a homeomorphism for theorem 1.13 to hold, it only needs to be onto, continuous and open, which leads to the concept of
Definition 1.15 (Semiconjugacy).
The maps and are said to be semiconjugate if there exists a map onto, continuous and open, such that .
That been stated, it can be shown that theorem 1.13 holds with “conjugate” replaced by “semiconjugate”.
Chapter 2 Regular-Reversal quadratic maps
In this chapter the work of [Solís and Jódar(2004)] will be discussed, as the main precedent to the present work, where an application of bifurcation theory and chaos in discrete dynamical systems is analyzed to determine the conditions in which chaos arises for a family of unimodal maps, namely quadratic maps. If the asymptotic behavior and parametric dependency of quadratic maps is understood, the results can be generalized to more complex maps, starting to set a precedent for cubic maps and polynomial maps of higher order. The only original contribution in this chapter is given by some of the examples in the corresponding section below.
The main objective of this chapter is to analyze, for a given unimodal quadratic map, under what circumstances does the property of period doubling hold and when it disappears, which means that chaos is no longer a possibility, so as to be able to control the appearance (or not) of chaotic behavior by controlling the value of the parameter.
2.1 Unimodal maps
Consider a one-parameter family of real discrete dynamical systems of the form with iteration function , , , intervals, and smooth on and .
Definition 2.1 (Unimodal Map).
An iteration function (map) , with an interval, is called unimodal if it is smooth111This requirement follows from the computability of the Schwarzian derivative. (at least ) and with a unique maximal point. By extension, we call a family of maps unimodal if every member of the family is unimodal.
In the following, let be unimodal and with a fixed point .
Definition 2.2 (Eigenvalue function).
The eigenvalue function corresponding to the fixed point of the one-parameter family of maps is
| (2.1) |
defined for the values of for which the fixed point exists.
That is, the eigenvalue function gives us the multiplier of a fixed point as a function of the parameter of the family of maps. The eigenvalue function for fixed is simply the derivative of the map at the fixed point, i.e. the “traditional” multiplier. Notice that is continuous since is smooth. The fixed point is asymptotically stable for the values of for which , agreeing with our notion of a hyperbolic attracting fixed point.
Definition 2.3 (Region of Type 1).
The region of type 1 of is the set .
The region of type 1 of is the set of values of the parameter for which is a stable fixed point, along with its corresponding values of the eigenvalue function (multipliers). Therefore, the region of type 1 of a fixed point is its stability region as a function of the parameter . Since is smooth on both and , regions of type 1 consist of the countable union of connected subsets of . In general, we can define regions of higher type respect to the fixed points of the -th iteration of , i.e. the periodic points of prime period .
Definition 2.4 (Region of Type ).
If the family of maps , not necessarily unimodal, has an isolated periodic point of prime period , , then we can define the -th eigenvalue function as
| (2.2) |
Notice that refers to the first derivative of the -th iteration of and not the -th derivative of , as in standard calculus notation. The set is denoted as the region of type for .
The region of type of the periodic point is simply the region of type 1 of the function , of which is a fixed point. Moreover, by the chain rule,
so that it suffices that the absolute value of the product of the multipliers along the orbit of a periodic point is less than one in order to guarantee the stability of . Nonetheless, the existence of periodic points of period is not assured even for a unimodal map in general. Thereafter, the periodic point, , when it exists, is stable in an interval of the parameter with ends given by and . We suppose here that the function exists, at least locally, i.e. for all in an interval around a value of interest associated with the periodic point , according to the Inverse Function Theorem [Protter(1998)]. A simpler way of saying the latter is that will be attracting as long as , but notice that depending on the form of the graph of , the eigenvalue function may never leave this region therefore undefining for . Another way of looking at the region of is that the type of a given region is given by the order of the attracting -cycle existing in that region.
Now, since a point can only belong to a unique and specific region of a given periodic point, we can state the following definition.
Definition 2.5 (Regular and Reversal Maps).
Let be an interval. A map is called a regular map (respectively, reversal) in if it has an isolated periodic point such that its associated eigenvalue function has the property that if , with , then the type of the region containing the point is higher (respectively, lower) or equal than the type of the region containing the point .
In a regular map, as the value of increases, new attracting periodic points of progressively higher periods are created and, therefore, the type of the regions associated with these points equally rises (or at least does not decrease) progressively. On the other hand, in a reversal map, as the value of the parameter increases, higher order attracting periodic points disappear and only progressively lower order attracting periodic points remain, so that the type of the associated regions decreases correspondingly.
Example 2.1 (The Logistic Map).
The logistic map is a classical example of a regular map, since it satisfies the definition 2.5 in the interval , with . Its eigenvalue function is
for the nonzero fixed point . In table 2.1 some type regions are shown for the logistic map.
| Type of the region () | Interval [] |
|---|---|
| 1 | (0, 3) |
| 2 | (3, 3.449489…) |
| 3 | (3.449489…, 3.544090…) |
| 4 | (3.544090…, 3.564407…) |
| 5 | (3.564407…, 3.568759…) |
| 6 | (3.568759…, 3.569692…) |
| 7 | (3.569692…, 3.569891…) |
Example 2.2 (The Exponential Map).
Consider the one-parameter family of maps , with . In figure 2.1 are shown three important cases of the graphs of the family. When the map has a repelling fixed point because of theorem 1.1 (see figure 2.1c); when , and , so that -1 is a non-hyperbolic fixed point and, by theorem 1.11, there is an interval around in which this fixed point exists and is unique; and when the fixed point is attracting also because of theorem 1.1. Therefore the fixed point of undergoes a change in stability from unstable to stable as the value of increases from to , with the bifurcation value being exactly . Also, it is easy to see that is an increasing function that is concave upward if and concave downward if (see figure 2.2). Thus has two fixed points at the points and shown in figure 2.2 when and they disappear when increases to . Since we know that has a single fixed point, these must in fact be periodic points of period 2. Note that, as these period 2 points lose their “attractiveness”, the fixed point acquires it, so that there is an “exchange of stability”. Also note that, at , , so that we can apply theorem 1.12 formally prove that a period doubling bifurcation effectively takes place. The bifurcation diagram of this example is shown in figure 2.3, where it is easy to see that it satisfies the definition of a reversal map.
Definition 2.6 (Regular-Reversal Map).
Let be a nonempty interval in . A map is called a regular-reversal map in if the interval can be decomposed in two nonempty disjoint intervals and such that the map is regular in and reversal in .
In a regular-reversal map, there is a value of the parameter , say , that divides the interval of the family of maps in two, in one of which it is regular and another one in which it is reversal. A regular-reversal map is interesting because it shows exactly the parametric range where periodicity takes place. So far, it has been difficult to determine for which maps this phenomenon occurs. It is for this reason that precedents have focused only in quadratic polynomial maps [Solís and Jódar(2004)], but also because quadratic polynomials give a complete family of unimodal maps with negative Schwarzian derivative and they represent the simplest example of unimodal maps. In this work, we will focus in higher order polynomial maps, particularly cubic, and we will be able to determine a simple way to create regular-reversal polynomial maps of any degrees. We can also define reversal-regular maps in a similar way.
2.2 Quadratic maps
Consider the family of quadratic discrete dynamical systems given by
| (2.3) |
where , and are smooth functions of the parameter and for all .
The fixed points of this system are obtained by solving the general quadratic equation
and they are given by the general quadratic formula
We will assume that in order to have two distinct real fixed points. We can simplify the system introducing the change of variable
where and . Note that is a fixed point of the system (2.3). We obtain the following modified discrete dynamical system [Solís and Jódar(2004)]
| (2.4) |
where . The map is unimodal. Without loss of generality, we will restrict to the case of positive and continuous for .
Under the above transformation, the fixed points of the new system are
| (2.5) |
so that the eigenvalue function for each of them is given by
The corresponding regions of type 1 are
There are also two periodic points of period two, which we can calculate by solving
We can factor from the expression since it corresponds to the fixed points (order 1) and we then solve for the roots of
and, again by the general quadratic formula, we have
with eigenvalue function
and, therefore, corresponding type 2 region given by
Moreover, there are periodic points of period with region of type given by
where the sequence is convergent, and converges to the limiting value [Solís and Jódar(2004)]. The sequence is related to the one described by [Feigenbaum(1978)] in [Feigenbaum(1978)] and converges exponentially. Table 2.2 summarizes this results.
| Type of the region () | Interval [] |
|---|---|
| 1 () | (-3, -1) |
| 1 () | (-1, 1) |
| 2 () | (1, |
| () |
The values of the original Feigenbaum sequence, , correspond to values of the points of intersection of the function with the the values of the sequence . Stated otherwise, the sequence results from evaluating the Feigenbaum sequence in the function . When the function is the identity function, we can recover the original Feigenbaum’s sequence.
Feigenbaum’s sequence, , describes the values of the parameter for which the period doubling bifurcation takes place in a unimodal map. The sequence converges exponentially. The values for a given iteration function –of the following form– are given implicitly by the definition of [Feigenbaum(1978)] in its original paper [Feigenbaum(1978)], for a map as
for the periodic point of period , , where we must solve for . In table 2.4 the values of a Feigenbaum sequence are shown for a specific iteration function.
A simple way to estimate the values of Feigenbaum sequences is starting from Feigenbaum’s constant, , and the first calculated values of the sequence , through the approximation
| (2.6) |
This estimate is used as a predictor for the next value of the sequence and tends to be more precise as .
2.3 Types of periodic points and chaos
By the known properties of the sequence and Sarkovskii’s theorem, we can state the following propositions [Solís and Jódar(2004)]:
Proposition 2.1.
If the function is bounded from above by , the map can only have periodic points of period with . Moreover, if the bound is given by , the system is not chaotic and it can only have periodic points of periods of powers of two.
It is worth remarking that there are unimodal maps that have no periodic points with periods greater than one, in which case the bifurcation diagram consists of a single branch given by the fixed point.
Proposition 2.2.
Suppose that the function is bounded. If the supreme of the function lies within the interval , then the map only has periodic points of period with and therefore, it is not chaotic.
In order for the system to be a regular reversal map for in some interval , it is necessary that the function has the following property: there exist in such that for some where is the -th element of the sequence . It is relatively easy to construct functions with the aforementioned property, which means there is a large class of unimodal maps that are not chaotic.
Beyond the value , there are other regions that divide the positive quadrant of . The upper region of , that is, the region , can be defined as the region of type and corresponds to divergent orbits. It is not clear how this region is subdivided and, in general, this is an open question.
2.4 Quadratic examples
The goal of the following examples is to show that regular reversal maps are a very broad family and that it is relatively simple to construct maps with specific –ad hoc– bifurcation diagrams with or without chaos.
Example 2.3.
Let be a quadratic unimodal function given by
with . Note that is the maximum value and is a root of (the other root is ). We are interested in knowing the periodic point structure of the map as we vary the parameter . Fixing , we then have the following cases:
-
•
If , we have a bifurcation diagram consisting of a single straight line coincident with the axis.
-
•
If , we obtain a regular reversal map with a bifurcation diagram consisting of a straight line, as the previous case, but with a closed asymmetrical loop homeomorphic to a circle (see figure 2.6a on page 2.6a). The lack of symmetry stems from the fact that the upper branch of the loop corresponds the the periodic point and the lower branch to . The diameter of the loop is calculated by finding the two intersections of with the line and it is equal to . The graph of is shown in figure 2.5a at page 2.5a.
- •
- •
- •
Example 2.4.
If we modify choosing it as periodic “unimodal”222Formally, this cannot be a unimodal function, since it is a requirement for a unimodal function to have a single maximum of the form
| (2.7) |
with ), the new bifurcation diagrams will be very similar to those of the last case, except that they will be periodic. In figure 2.7 at page 2.7, the graphs of are shown, and in figure 2.8, are the corresponding bifurcation diagrams, both for the corresponding parameter cases of figures 2.5 and 2.6, respectively.
Example 2.5.
If we introduce an asymmetry factor to the function 2.7 we can achieve an asymmetrical bifurcation diagram with the different lobes having different amplitudes. One way to accomplish the above is defining as
| (2.8) |
Varying the parameter we obtain graphs of like the ones shown in figure 2.9 at page 2.9, with its respective bifurcation diagrams, shown in figure 2.10 at page 2.10. For small values of we have two lobes (or “humps”) in the function , with its local maxima under the value , which translates into a bifurcation diagram with four asymmetrical lobes and of different amplitude (see figure 2.10a at page 2.10a). However, due to asymmetry, rising the value of to will only produce the appearance of nested loops in two of the right lobes of the bifurcation diagram (see figure 2.10b), which repeats itself as we continue to raise the value of (see figure 2.10c).
Example 2.6.
A similar effect to the one of the last example can be achieved by taking as a polynomial of fourth grade in , with its coefficients chosen in such way that it is not unimodal. For example,
| (2.9) |
with . The function has three local extremes (see figure 2.11 at page 2.11). In this case, the bifurcation diagram is shown for different values of the parameter in figure 2.12 at page 2.12 where, for progressively higher values of we progressively obtain lobes that stem from the branches corresponding to the periodic points. In this case the diagram is no symmetrical nor periodic because the function is neither.
With the above examples we have shown we can design discrete dynamical systems in order to have specific bifurcation diagrams “on demand” with desired properties. The fundamental fact that allows us to do the latter is to note where the function intersects the values of the sequence . The detailed form of is irrelevant as long as it does not cross a value of the sequence, for it only affects the specific form of the “branches” of the periodic points in the bifurcation diagram, but it does not alter the fundamental structure of the set of periodic points.
Chapter 3 Quadratic maps revisited
Now, we will restate and expand the work of [Solís and Jódar(2004)] in a way that will allow us to gain greater insight into the way we can control the appearance of chaos in a parametric family of quadratic maps and we will be able to generalize these results in a natural way for cubic maps and, later on, for general polynomial maps of -th degree.
3.1 General quadratic form
Consider the quadratic iteration function as stated by the
Definition 3.1 (General Quadratic Map).
The General Quadratic Map (GQM) is defined by
| (3.1) |
where
| (3.2) |
and , and are real coefficients that depend upon the parameter and are . is called the Fixed Points Polynomial (of the general quadratic map).
It is clear that any quadratic map can be put in this form. Although the minus sign in the definition of may now seem “unnatural”, its function will be clear shortly111Check page 3.2 below..
Evidently, the zeros of are the fixed points of –thence the name of – since, if is a root of , then implies .
In the following, for simplicity, we will drop the subscript “2” from to denote the CQM. Since the derivatives of are
| (3.3) | ||||
straightforward calculations lead us to
Lemma 3.1.
has the following two roots and critical point, respectively,
| (3.4) | ||||
with corresponding extreme value
And, finally, it has the fixed points
| (3.5) | ||||
3.2 Linear factors form
Suppose we know both roots of and they are , which exist by the fundamental theorem of algebra. Then, by the factor theorem, we know that is a factor of and we can rewrite it as
| (3.6) |
Therefore, we can also rewrite as
| (3.7) |
Since , then we can define , where
| (3.8) |
is the sign function and ; this definition will be useful for our purposes further below.
With the above, we have
Definition 3.2 (Linear Factors Form of the Quadratic Map).
The GQM can be rewritten as
| (3.9) |
where and are smooth functions of the parameter . We call the Linear Factors Form of the Quadratic Map (LFFQM).
Again, straightforward calculations lead us to
Lemma 3.2.
(3.9) has first and second derivatives with respect to given by
| (3.10) | ||||
The roots of are
| (3.11) | ||||
and the only critical point is at
with corresponding extreme value
Here the apparently “unnatural” minus sign of the definition 3.1 of the fixed points polynomial will show its usefulness: for purely aesthetic reasons, in the definition of the LFFQM, we want that whenever and which, without loss of generality, is accomplished by arbitrarily introducing the mentioned sign.
3.3 Canonical form
We will introduce a canonical form of the quadratic map by requiring both fixed points in the linear factors form (3.2) to be real or, equivalently, in the general form (3.1); also, we will require an “amplitude”, , equal to unity222Recall is the function preceding the product of the linear factors of the fixed points in the LFFQM defined in Eq. (3.2)., a fixed point mapped to zero333This is the reason we require both fixed points to be real, otherwise we would need a complex transformation to map the complex fixed point to zero, which is beyond the scope of this work. and the other fixed point being mapped to
| (3.12) |
The latter is accomplished by taking one of the following linear transformations
| (3.13) | ||||
Without loss of generality, we will take transformation and call it simply . Since is evidently linear, its inverse can be calculated directly as
| (3.14) |
Performing the corresponding calculations on the LFFQM we obtain
Definition 3.3 (Canonical Quadratic Map).
The map defined by
| (3.15) |
where is a parameter, is called the Canonical Quadratic Map (CQM).
We have thus mapped the real fixed points and . This form is much simpler than the previous ones and it has the advantage of putting the whole parametric dependence onto the single non-constant fixed point through (3.12). Implicit in this dependence is the parametric dependence of the corresponding original roots of the LFFQM (3.9), so we can explicitly give this dependence if desired. Moreover, the fact that is actually a topological conjugacy between and will be proved in chapter 5, by which we can state that the study of stability and chaos in quadratic maps with real fixed points can be summarized by the study of the Canonical Quadratic Map, which yields its importance.
With some straightforward calculations, we easily arrive at the following
Lemma 3.3.
The Canonical Quadratic Map has the roots
| (3.16) | ||||
and has derivatives
| (3.17) | ||||
so that its critical point is at
with value
| (3.18) |
We see then that the second root of the canonical map, , is always one unit to the right of the non-constant fixed point . Also, the critical point is always a maximum and is located in the middle point between 1 and , with its value being simply the square of the point.
3.4 Stability and chaos in the canonical quadratic map
3.4.1 Fixed points
The stability of the fixed points is given by the absolute value of its multiplier, which must be less than one. Making use of the defined eigenvalue function, we have from (3.17) that for the fixed point ,
| (3.19) |
so that it is straightforward to see that if
| (SC0) |
then is a stable fixed point. We will call inequality SC0 the stability condition for the fixed point . We see then that the stability of actually depends on the value of . In terms of the original fixed points of the quadratic map, and , what counts is the separation between them scaled by the factor , as can be seen from equation (3.12).
On the other hand, for ,
| (3.20) |
so that is stable as long as
| (SC1) |
We call the value our first bifurcation value. Notice that is also a bifurcation value since it is the limit of the stability condition for , given by (SC0). We summarize this results in the following
Proposition 3.4.
Let and be the two fixed points of the canonical quadratic map as defined above. Then
-
•
if , zero will be an asymptotically stable fixed point;
-
•
if , will be an asymptotically stable fixed point.
3.4.2 Periodic points of period two
To determine the period two periodic points we must of course solve
Some basic algebra shows that
Therefore, we must solve
Clearly, we must factor out and since we know those are the fixed points of the original map, and we are looking (only) for period two points. Performing the factorization, we get
| (3.21) |
Solving the latter equation we get to the following
Lemma 3.5.
The canonical quadratic map has a unique 2-cycle given by
where the superscript will indicate the period of the periodic point.
Proof.
It is straightforward to see that and , therefore proving that both and are period 2 periodic points of and the orbit is indeed a 2-cycle. The uniqueness comes from the fact that both are algebraic roots of the second degree equation (3.21) and every such period 2 periodic point must be so. ∎
To analyze the stability of this 2-cycle, we must calculate its own eigenvalue function. The derivative of is
so that evaluating in and we have
| (3.22) |
so that, not surprisingly, both points of the 2-cycle have the same stability criteria. The stability region is then determined by
whose solutions lead us to
Lemma 3.6.
The unique 2-cycle of the canonical quadratic map is asymptotically stable whenever one of the following inequalities is met
where is the nonzero fixed point of the CQM.
So, we see that the stability region is divided in two in terms of the value of the original fixed point . As a reference, . We call this value the second bifurcation value of the parameter for the canonical quadratic map. Notice again that is also a bifurcation value.
It is worth noting that in the case of the first inequality of lemma 3.6, it is the zero fixed point that produces a bifurcation, while for the second one, it is that bifurcates to produce the 2-cycle. The latter can be seen from lemma 3.5 since, for , the 2-cycle is not real and it only begins to be so at the bifurcation values ; when , both and are zero (and equal to the zero fixed point) and thereon separate as the two branches of the square root function (positive and negative). On the other hand, when , both and are also 2 in value and thereon separate in the same way as before.
We can now prove that, indeed, what happens at the bifurcation values of are period doubling bifurcations.
Proposition 3.7.
Let444The dot represents derivative with respect to the parameter to avoid confusion with the apostrophe representing derivative with respect to . (respectively, ) and . Then the zero fixed point (respectively, the fixed point) of the CQM undergoes a period doubling bifurcation precisely when .
Proof.
According to theorem 1.12 we must prove that:
-
1.
for all in an interval around . This is trivially true since the zero fixed point does not vary with the parameter.
-
2.
. Which follows directly from eq. 3.19.
-
3.
. From (3.22) we see that this is true since by hypothesis.
The case for the period doubling bifurcation of the other fixed point when is analogous. ∎
The requirement will be understood more clearly in the examples below. Likewise, we can prove that, in turn, the fixed point (or ) of —periodic point of period 2 for — also undergoes a period doubling bifurcation associated with the second bifurcation value when .
Proposition 3.8.
Let and . Then the fixed points and of the undergo a period doubling bifurcation precisely when .
Proof.
Consider first the case for given by lemma 3.5. According to theorem 1.12 we must prove that:
-
1.
. Which follows directly from eq. 3.22 since by hypothesis .
-
2.
for all in an interval around . Theorem 1.11 allows us to affirm this since by the above point .
-
3.
. Since, by the chain rule,
using eq. (3.22) and the fact that is a fixed point of . The partial derivative with respect to is then
since and by hypothesis.
The cases for and are analogous. ∎
3.4.3 Periodic points of higher period
Some more algebra shows us that
| (3.23) |
with derivative
| (3.24) |
So that, in order to find the periodic points of period three we must, in principle, explicitly solve
| (3.25) |
which is an seventh degree polynomial in . Even factoring out the other known fixed point of , , we are still left with a sixth degree polynomial on which is, in general, impossible to solve explicitly by algebraic means. Nevertheless, for specific polynomials we can determine the roots numerically.
Much less can be said about determining higher period periodic points. However, we may still approximate the stability regions of these -cycles by numerically determining the bifurcation value of the parameter , which in this case we have done by plotting very high-precision bifurcation diagrams (see Appendix A). One way to determine an approximation to the value of the infinite period bifurcation value, after which the onset of chaos takes place, is through equation (2.6), which yields a very rough estimate of . Some bifurcation values for the CQM found in this work are shown in table 3.1; moreover, table 3.2 shows the the stability conditions for each periodic point,up to period ; the boundary of the stability bands correspond to the bifurcation values.
| 0 | 0 |
| 1 | 2 |
| 2 | |
| 3 | |
| 4 | |
| 5 | |
| 6 | |
| 7 | |
| Period | Periodic Points | Stability Condition |
| 1 | ||
| 1 | ||
| 2 | ||
| 4 | ||
| 8 | ||
| 16 | ||
| 32 | ||
| 64 | ||
| ⋮ | ⋮ | ⋮ |
| - |
3.4.4 Fixed points with multiplicity
In the CQM, the only way we can have multiplicity in the fixed points is when . In this case, takes the form
| (3.26) |
Its derivative is then
| (3.27) |
and, therefore, , so that it is a nonhyperbolic fixed point, according to definition 1.6.
Proposition 3.9.
Let be the CQM with a single fixed point with multiplicity of two, as in eq. (3.26). Then the single fixed point is semistable from the right.
Proof.
It is straightforward to see that and then for all . Therefore, in particular, and then , so that, by the latter, we cannot use theorem 1.7. However, theorem 1.9 tells us that the fixed point is semistable from the right. Indeed, as is concave downward at , then is decreasing in a small neighborhood about zero and, since and is continuous, we have that there exists , so that for and for . Consequently, is unstable from the left and semistable from the right, according to definition 1.9 and the proof of theorem 1.1 [Elaydi(2000)]. ∎
It is worth noting that, although the case of multiplicity in the fixed point of the CQM is possible, the case of complex fixed points is not, since the CQM requires one fixed point to be zero, then its complex conjugate —which should also be a fixed point— is itself.
3.4.5 Regular-reversality and chaos in the canonical quadratic map
Let us recall, that the definition of transforms, through , the parametric dependence of the original fixed points of the general quadratic map into a single parameter-dependent function . Recall then that,
Therefore, assuming a constant value of , we see that the value of is exactly the difference between the two original fixed points, and , of the General Quadratic Map (which we forced to be real). We can then interpret the found values of the sequence as determining the
Definition 3.4 (Stability Bands of a Quadratic Map).
Graphically, we see then that stability bands run along the values of and as functions of the parameter and, as long as the value of one fixed point does not cross a the limit of a band of the other fixed point, there will be no changes in the periodic points structure of the map, i.e. there will only be bifurcations when the value of one fixed point crosses the limit of a stability band of the other fixed point. If were not held constant, it would simply act as a “modulator” for the width of the bands along the parameter .
Stability bands for and in the CQM are defined analogously.
Definition 3.5 (Stability Bands of the CQM).
Notice that that the part of the -th stability band in the upper (respectively, lower) semi-plane corresponds to the stability condition of period attracting periodic points associated with period doubling bifurcation of the (respectively, zero) fixed point.
Finally, we must remark that the width of these stability bands —given by the values of the sequence — is not equal for polynomial maps of different degrees, as we shall see in the next chapters.
By the above, the big picture is that depending on the absolute value of the one-parameter family of canonical quadratic maps can present or not attracting periodic points of different periods or even become chaotic. In particular, as long as we will only have attracting periodic points of period less or equal to . We can thus restate the propositions 2.1 and 2.2 given in section 2.3 of chapter 2, originally by [Solís and Jódar(2004)] ([Solís and Jódar(2004)]).
Proposition 3.10.
If the absolute value of the fixed point is bounded from above by , the GQM map is not chaotic and can only have periodic points of period with . Moreover, if the bound is given by , the system is still not chaotic and it can only have periodic points of periods of powers of two.
And more specifically,
Proposition 3.11.
Suppose that the absolute value of the fixed point is bounded. If the supreme of the function lies within the interval , then the map only has periodic points of period with and therefore, it is not chaotic.
In order for the system to be a regular reversal map for in some interval , it is necessary that the continuous function has the following property: there exist in such that for some where is the -th element of the sequence of bifurcation values . It is easy to construct functions with the aforementioned property, which means that constructing regular-reversal canonical quadratic maps is also easy. Furthermore, knowing beforehand the properties of a desired bifurcation diagram, i.e. the fixed points, the values for desired bifurcations, etc. we can construct functions that accomplish the desired properties in the bifurcation diagram.
3.5 Quadratic examples with the canonical map
Recall that the CQM is
where
| (3.30) |
Therefore, we can “parametrize” the map by defining either directly or going back to the parametric dependence of the original LFFQM from which was defined using the transformation , i.e. defining , and . We will look at both interpretations at the same time in the following examples.
3.5.1 Linear fixed points
We will begin with the simplest case: linearly varying and and . Consider,
The graphs of these parametrizations of the fixed points are shown in figure 3.1. Clearly, and its graph is also shown in figure 3.1.
The bifurcation diagram for this specific parametrization of the fixed points of the CQM is shown in figure 3.2. In this bifurcation diagram we can see that it is precisely at the values of for which crosses limits of the stability conditions of (or, equivalently, where crosses the bifurcation values) that bifurcations take place. Notice that, since , it is the fixed point that undergoes period doubling bifurcation.
In fact, it does not matter what the actual values of the slopes of and are: as long as their difference is the same (has the same slope) they will produce exactly the same bifurcation diagram, since the parametric dependence of will be unchanged (see eq. 3.30). Making precise calculations of the bifurcations diagram near -values of interest we can determine numerically the rest (or at least some more) of the bifurcation values forming the sequence , but this is best accomplished using the parametrization of example 3.5.5 below.
3.5.2 Quadratic fixed point
Redefining the parametric dependencies of and we can achieve regular-reversal maps as described in definition 2.6 in page 2.6. For example, leaving in the above example 3.5.1, but redefining with a quadratic parametric dependence on as
| (3.31) | ||||
We immediately see that zero is a root of the parametric form of , and the other root is . Since it is exactly like the defined Linear Factors Form defined for the quadratic map in subsection 3.2 on page 3.2, we know that its maximum lies at
with corresponding value
Using, for example, the values, and we know then that and, therefore, (this functions are shown in figure 3.3a. As for it is then
where the maximum is clearly located at with a value of . The graph of is shown in figure 3.1b. Also, crosses the first bifurcation value () two times at and , and it is easy to see that for those values of , so that by proposition 3.7, we expect period doubling bifurcations at this values. Also, since , according to proposition 3.11 we know that we will only get one such bifurcation leading to stable 2-cycles in the interval , and then going back to a unique stable fixed point (i.e. stays in the 1-stability band in the interval ). Thus, this example complies with our definition of a regular-reversal map. The latter predictions are confirmed in the bifurcation diagram shown in figure 3.4.
Now, increasing the maximum value of we can accomplish more bifurcations in the diagram. This can be achieved by increasing the slope of to , since for this value . The corresponding graph of is shown in figure 3.5a. And, also, the respective bifurcation diagram is shown in figure 3.5b, where we can see that the system undergoes a cascade of period doubling bifurcations in the central part. We thus produce another regular-reversal map with more nested regions of different types.
If we continue increasing the maximum value of we can achieve a chaotic region within the regular-reversal map. Any value above the used last above will suffice. Say we take . The graph of for this case is shown in figure 3.6 and we then get the bifurcation diagram from figure 3.7, which displays a chaotic region in the central part, as stated. This is still a regular-reversal map.
3.5.3 Sinusoidal fixed point
We can construct a periodic regular-reversal map by making one of the fixed points have a sinusoidal parametric dependence on . That is, we leave as before and define
| (3.32) |
where we can tweak the sinusoidal form of by means of the constant which is just an amplitude (and translation, at the same time). Choosing we can get a periodic regular-reversal map with chaotic regions in the periodic central regions. The graphs of the roots are shown in figure 3.8 and the bifurcation diagram is shown in figure 3.9.
3.5.4 When zero loses its stability
It is important to point out that if becomes negative, i.e. the difference between and becomes negative (assuming again ), then we enter the stability region of the fixed point at zero, , in the canonical quadratic map (and in the general quadratic map). Thus decreasing values of will cause bifurcations to appear this time from the zero fixed point, as the stability bands are crossed by (correspondingly, in the general quadratic map). We can achieve the latter by defining again a linear as
with, say and . Then we get to cross through all the stability bands of in the interval of the parameter . This is shown graphically in figure 3.10 along with its corresponding graph of .
The bifurcation diagram in this case starts, from left to right, being chaotic, then reversal, then switching to stability of (correspondingly, ) and then being regular until chaos is reached again, as crosses the stability bands. This is shown in figure 3.11.
3.5.5 Asymptotic fixed point
Since we know that the onset of chaos occurs at (as calculated through approximation formula 2.6), we can force the moving fixed point of the CQM to asymptotically approach this value in order to visualize more clearly the bifurcation process that takes place very near this value. Now, for symmetry reasons, it will be best if we make the fixed point to be the one that bifurcates, instead of ; we do this by making (remember is the interval where the period doubling bifurcations cascade occurs for ). So, we define
| (3.33) |
where and . We choose for this example, for the above reasons, . The corresponding bifurcation diagram for this case is shown in figure 3.12, upper panel, where we can clearly see the bifurcations take place up to the attracting periodic point of period (64-cycle); the lower panel in the same figure shows the value of the fixed points, which shows the asymptotically varying . The scale of the parameter value is the same in both panels to allow for comparison; notice that, when the fixed point curve crosses the marked values, period doubling bifurcation takes place in the bifurcation diagram, as expected. By making precise calculations based on this diagrams we could determine the bifurcation values up to an arbitrary -cycle, as we did for up to in this work.
Chapter 4 Cubic maps
As stated before, in chapter 3 we furthered the results of chapter 2, mainly based on the work of [Solís and Jódar(2004)], and restated them in the frame of a new formulation that would allow for generalization. Well, the first such generalization comes precisely for third degree polynomial maps, i.e. cubic polynomials.
4.1 General Cubic Map
Consider a cubic map (iteration function) in its general form, as stated by the
Definition 4.1 (General Cubic Map).
The General Cubic Map is defined by
| (4.1) |
where
| (4.2) |
is called the Fixed Points Polynomial of . All the coefficients and are functions of the parameter .
It is evident that any cubic map can be put in this form by adjusting the corresponding values of the coefficients in the fixed points polynomial. By the fundamental theorem of algebra, we know that (4.2) has three roots, by which the GCM has three fixed points. The roots of are then [Spiegel et al.(2009)Spiegel, Lipschutz, and Liu]
| (4.3) | ||||
where , , , and
| (4.4) | ||||||
The coefficients of eq. (4.2) and its fixed points are related by
| (4.5) |
is called the discriminant and we have three cases:
-
•
If then one fixed point is real and the other two are complex conjugates.
-
•
If then the three fixed points are real with at least two of them equal.
-
•
If then all fixed points are real and distinct.
It is the last two cases (real fixed points) that will interest us most for the time being. Suppose in particular that . Then, we can write [Spiegel et al.(2009)Spiegel, Lipschutz, and Liu]
| (4.6) | ||||
where .
4.2 Linear Factors Product Form
Since we know that (4.2) has three roots, by the factor theorem we can rewrite as
| (4.7) | ||||
where , ; and then rewrite (4.1) as
Definition 4.2 (Linear Factors Form of the Cubic Map).
Let be a general cubic map with three fixed points, , and . We can write as
| (4.8) | ||||
where all and are functions of the parameter . We call the Linear Factors Form of the cubic map (LFFCM).
The LFFCM has the derivatives
| (4.9) | ||||
By equating the first derivative to zero we can obtain the two critical points
| (4.10) | ||||
and, evaluating these critical points in we have
| (4.11) | ||||
so that necessarily one critical point is a maximum and the other one is a minimum, unless they are both equal, in which case we have a saddle point.
4.3 Canonical Cubic Map
As in the case for the Canonical Quadratic Map, we will apply a linear transformation to (4.8) so that one fixed point is mapped to zero and the “amplitude” of the linear factors term is unity. We know that, since is cubic, at least one fixed point is real, so we can map this fixed point to zero. This can be accomplished by one of the following transformations
| (4.12) |
by taking if is real. Without loss of generality, we will assume is real and apply with the plus sign calling it simply , so that, once the calculations are done, we get
Definition 4.3 (Canonical Cubic Map).
The Canonical Cubic Map (CCM) is defined by
| (4.13) |
where it has been stressed out that both fixed points and depend upon the parameter .
So that, if then
| (4.14) |
and if
| (4.15) |
From the easily verifiable calculations, it can be shown that the relation between the roots of the linear factors form of the cubic map and the canonical cubic map is given by
Lemma 4.1.
The fixed points of the linear factors form of the cubic map and the canonical cubic map are related by
| (4.16) |
We have then again reduced the parametric dependence to only two functions of the parameter : and . Notice is a homeomorphism between the domains of both maps; this will help us in chapter 5 to prove that the Linear Factors Form and the Canonical Form of polynomial maps are actually topologically conjugate, which in turn means that the stability and chaos properties are preserved between the maps, which allows us to determine stability properties for any cubic map by analyzing only the CCM.
Now, one first immediate result is
Proposition 4.2.
The Canonical Cubic Map has three roots:
| (4.17) |
Moreover, if
-
•
If and , we will have a single root different from zero with multiplicity of two, and equal to , i.e. the midpoint between the two non-zero fixed points.
-
•
If and , zero will be the only root, with multiplicity of three.
-
•
If , we will have two real and distinct non-zero roots.
-
•
If , zero will be the only real root of the CCM and the other two roots will be complex conjugates.
And if the three roots are always real and distinct, except when when has multiplicity of two.
As for the derivatives of the CCM, they are
| (4.18) | ||||
from which we can determine that the critical points are
| (4.19) | ||||
So that, evaluating at the critical values of and we have thatwe have that
and, then, we have as expected one maximum and one minimum unless both are equal, in which case we have a saddle point. From this we have proved
Proposition 4.3.
Consider the CCM as defined above. Then, if
-
•
If , both critical points are real and distinct, one corresponding to a maximum and the other to a minimum.
-
•
If , both critical points are equal and correspond to a saddle point.
-
•
If , the two critical points are complex conjugates, so we do not have any real critical points.
And if ,
-
•
If , both critical points are real and distinct, one corresponding to a maximum and the other to a minimum.
-
•
If , both critical points are equal and correspond to a saddle point.
-
•
If , the two critical points are complex conjugates, so we do not have any (real) critical points.
4.4 Stability and Chaos in the Canonical Cubic Map
Next, let us determine the stability of the periodic points of the CCM. As we will see further below in chapter 5, this analysis will suffice for any cubic map, by means of topological conjugacy. However, we can only explicitly give this for the fixed points.
4.4.1 Fixed Points
We already know, by construction, that the fixed points of the CCM are and . While the first is constant always, the other two fixed points are set to be functions of the parameter . By evaluating in we get the eigenvalue functions. For we have
So that, the stability condition for this fixed point is
| (4.20) |
We can draw some conclusions from this. In order for zero to be a stable (attracting) fixed point one must have:
Lemma 4.4.
The following are sufficient conditions for the asymptotic stability of the zero fixed point of the Canonical Cubic Map:
-
•
in magnitude, ; and
-
•
if , and must have different signs; or
-
•
if , and must have the same sign.
Notice that the stability condition (4.20) states that the product of the distances from the other two fixed points to the zero fixed point must within the range , for positive [or for negative ], for the zero fixed point to be asymptotically stable.
The case of , , is not included in the discussion here since this would represent repeated fixed points (multiplicity), which will be discussed in section 4.6 below; likewise, in the remainder of this section we will avoid dealing with multiplicity of the fixed points.
Now, for , its eigenvalue function is
so that the stability condition for this fixed point is
| (4.21) |
This gives us the following
Lemma 4.5.
The following are sufficient conditions for the asymptotic stability of the fixed point:
If then
-
•
and must have the same sign; and
-
•
.
On the other hand, if ,
-
•
; and
-
•
if , then ; or
-
•
if , then ; or
-
•
if , then .
Again, notice that the stability condition (4.21) for can be translated as that the product of the distances between the other two fixed points and must be within the range for positive [or for negative ]. Also notice that when the bound may be negative even if or positive even if , therefore the usefulness of the distinction.
Finally, for we have the eigenvalue function
which gives the stability condition
| (4.22) |
which in turn produces the following
Lemma 4.6.
The following are sufficient conditions for the asymptotic stability of the fixed point of the Canonical Cubic Map:
If then
-
•
and must have the same sign; and
-
•
.
On the other hand, if ,
-
•
; and
-
•
if , then ; or
-
•
if , then ; or
-
•
if , then .
Likewise, we see that the stability condition (4.22) involves the product of the distances from the other two fixed points to the one whose stability is of interest. We will later generalize these “stability conditions” to functions of the parameter which are different for each fixed point, but of whose value depend the stability of not only the fixed points, but higher period periodic points also, through period doubling bifurcations. From the stability conditions for the three fixed points we have proved the following
Corollary.
A cubic polynomial map with three different real roots can only have a single attracting fixed point.
Proof.
Compare the stability conditions for the three fixed points. ∎
And,
Theorem 4.7.
If , then sufficient conditions for the stability of a fixed point of the canonical cubic map are
-
•
The product of the distances between each unstable fixed point and the stable one must be negative, which means one distance is positive and the other negative, which leads us to
-
•
the fixed point that lies between the other two will be stable, while the outer fixed points will be unstable, as long as
-
•
the product of the distances between each unstable fixed point and the stable one must be greater than -2.
4.5 Higher period periodic points
Although, as previously stated, in general, we cannot calculate the values of the periodic points of period 2 or higher, we can calculate for which values of the stability conditions above the fixed points undergo period doubling bifurcations. We will see in chapter 5 that these stability conditions can actually be generalized to something called the “Product Distance Function”, which depends on the parameter and is different for each fixed point. Example 4.4 below allowed us to determine the bifurcation values, , of the fixed points of the CCM up to some precision. The values obtained are shown in table 4.1. When the stability conditions of each fixed point crosses these values, bifurcations take place.
| 0 | 0 |
| 1 | 2 |
| 2 | |
| 3 | |
| 4 | |
| 5 | |
From these values, we can construct the analogue of the stability bands of the CQM for the CCM.
Definition 4.4 (Stability Bands of the CCM).
Notice, however, that in contrast with the stability bands of the CQM, the stability bands of the CCM cannot be plotted along the fixed points plots, at least not directly as just defined, but rather they must be represented in a separate plot for the stability conditions, as we shall see in the examples of section 4.7 further below.
4.6 Multiplicity of the fixed points
When multiplicity of the fixed points takes place in the CCM, without loss of generality, can take the following forms
| (4.24) |
with corresponding derivatives
| (4.25) |
and therefore, , for all three cases, so that we deal with nonhyperbolic fixed points, according to definition 1.6.
Proposition 4.8.
The stability of the fixed points of the CCM when they present multiplicity is, without loss of generality, as follows
-
1.
If , the zero fixed point is an unstable fixed point with multiplicity of two; and
-
•
if and
-
–
if it is semiasymptotically stable from the right;
-
–
if it is semiasymptotically stable from the left;
-
–
-
•
or if and
-
–
if it is semiasymptotically stable from the left;
-
–
if it is semiasymptotically stable from the right.
-
–
-
•
-
2.
If , this fixed point has multiplicity of two and it is unstable; moreover
-
•
if and
-
–
if it is semiasymptotically stable from the left;
-
–
if it is semiasymptotically stable from the right;
-
–
-
•
or if and
-
–
if it is semiasymptotically stable from the right;
-
–
if it is semiasymptotically stable from the left.
-
–
-
•
-
3.
If , the zero fixed point has multiplicity of three and
-
•
if , it is unstable;
-
•
if , it is asymptotically stable.
-
•
Proof.
Notice that
| (4.26) |
and
| (4.27) |
for all cases. Therefore
- 1.
- 2.
-
3.
If , the zero fixed point has multiplicity of three and we have that , and ; therefore, by theorem 1.7, if the zero fixed point is unstable and if it is asymptotically stable.
∎
4.7 Cubic examples
Here we will deal with specific parametrizations for the fixed points and in order to clarify the above findings and to demonstrate how we can construct bifurcation diagrams with specific predetermined properties with cubic maps. We will consider unless otherwise stated explicitly.
Example 4.1.
First, consider linear parametrizations for both and as
| (4.28) |
The result is plotted in the lower panel of figure 4.1, where we see that the middle fixed point is always, so we expect this to be the only stable fixed point, until the separation between this and the other points breaks the stability condition and the period doubling bifurcation cascade sets on. The corresponding stability conditions are shown in the middle panel of figure 4.1, where we confirm that the curve for is the only one within the stability band until where the curve crosses the barrier of -2 and drops from then, causing to bifurcate. The corresponding bifurcation diagram is shown in the upper panel of figure 4.1, where we confirm the stated above.
Example 4.2.
Next, based on the fact that the stability conditions depend upon the product of the distances between the stable fixed point and the other two fixed points, we take square root parametrizations of the fixed points and so that the product is linear. We take the parametrizations as
| (4.29) |
The plot of the fixed points is shown in lower panel of figure 4.2, where we see again that the middle fixed point is , so this will be the only attracting fixed point as long as its corresponding stability curve is within the stability region. In the middle panel of figure 4.2 we see that such stability curve goes out of said region at and continues to bifurcate (period doubling) thereof. The corresponding bifurcation diagram is shown in the upper panel of figure 4.2, which is very much like the one of example 4.1 above, but looks “stretched” since the progression towards bifurcation is now linear, instead of quadratic, and given by the stability conditions curves.
Example 4.3.
Now we will explore the full range of stability regions by making a linearly varying fixed point pass through the regions defined by the constant and a constant . We define then
| (4.30) |
We then obtain the plot of the lower panel of figure 4.3, for the selected range of interest of the parameter . In the middle panel of the same figure we can see the stability curves for the fixed points, where we see that initially, from left to right, all fixed points are unstable then, progressively, , and become stable, the latter one losing stability for still greater values of . The corresponding bifurcation diagram is shown in the upper panel, where we can see how first, the stable fixed point is , since it is the middle one, but begins in the chaotic region and goes “reversal” towards being stable; then, as crosses through zero, it becomes the middle stable fixed point and, when it in turn crosses the constant , this latter one becomes the stable fixed point, again loosing stability when crosses the stability band for .
Example 4.4.
This example is analogous to example 3.5.5 of chapter 3, where we had the fixed point of the CQM tend asymptotically to a value of approximately . To achieve the same qualitative phenomenon in the CCM, we can define the roots as
| (4.31) |
so that and the fixed points are symmetric about the zero fixed point; their product, which gives the stability condition for the zero fixed point, will then be which will asymptotically tend to , as ; this allowed us to calculate the bifurcation values of table 4.1. The lower panel of figure 4.4 shows the graphs of the fixed points, the middle panel shows the stability conditions for each fixed point and the upper panel shows the corresponding bifurcation diagram. This example allows us to determine the bifurcation values of the CCM with more precision; the values obtained are shown in table 4.1.
Chapter 5 n-th degree polynomial maps
Consider again a one-dimensional discrete dynamical system defined by
| (5.1) |
Where is a polynomial in one real variable with real fixed points and whose coefficients depend smoothly on the real parameter . Depending on the form of we have arrived previously to the definitions of the General, Linear Factors and Canonical Forms of the polynomial quadratic and cubic maps. Next, we will define precisely what we understand by the General and Canonical Maps of -th degree.
Definition 5.1 (General Polynomial Map).
The General Polynomial Map of -th degree (GPM-) is defined by
| (5.2) |
where
| (5.3) |
is called the Fixed Points Polynomial associated to .
It is obvious that any -th degree real polynomial in one variable can be put into the General Form by means of adjusting the value of the coefficient properly in the fixed points polynomial. This is the broadest class of real polynomials of finite degree.
It is clear that the roots of are the fixed points of . In the case of odd, the fundamental theorem of algebra guarantees the existence of at least one real fixed point. Let , be the roots of , then is a factor of by the factor theorem, therefore we can rewrite as
| (5.4) |
and then define the
Definition 5.2 (Linear Factors Form).
We can directly verify that for we obtain the corresponding Linear Factors Form of the Quadratic and Cubic Maps, described on pages 3.2 and 4.2, respectively. Once we know the fixed points of a map , it is straightforward to write its linear factors form. The motivation behind the factor is that we want that, for purely aesthetic reasons, if , the fixed points are real and , we have that , which is thus accomplished.
We will now restrict this set of polynomials to those whose fixed points polynomials have only real roots, i.e. maps with real fixed points only, though not necessarily distinct. Let us make this precise by the following
Definition 5.3 (Canonical Polynomials Set).
The Canonical Polynomials Set, denoted , is
| (5.6) |
where is the set of polynomials with real coefficients on the variable . Likewise, denotes where is the set of polynomials of degree with real coefficients on the variable .
The set has been our main work ground for the analysis in this work and, as it turns out, its elements can be put in a much nicer form, easier to understand.
We can further reduce the complexity of this set of maps by means of the transformation
| (5.7) |
where is a real fixed point of the map in its linear product form. Notice that is linear, therefore it has an inverse
| (5.8) |
is in fact a homeomorphism. We will drop the subscript when referring to the transformation for no specific degree. Applying this transformation to we reach the
Definition 5.4 (Canonical Polynomial Map).
The Canonical Polynomial Map of -th degree (CPM-) is
| (5.9) |
where
and are the fixed points of the corresponding linear factors form map of -th degree, (at least) being real.
It is clear from the definition that always. Notice also that the result from evaluating in the corresponding . We can easily prove that not only does the canonical map result from applying , but that the canonical map is in fact -conjugate to the linear factors form.
Proposition 5.1.
Let and and as defined above, having at least one real fixed point; let be this real fixed point, without loss of generality. Then is -conjugate to .
Proof.
It is clear that is a homeomorphism since it is linear. Then, we must only prove that , i.e. . We then have
| (5.10) | ||||
Where we have used and . ∎
This turns out to be very useful, since we know that topological conjugacy is an equivalence relation that preserves the property of chaos (see theorem 1.13 in page 1.13). This means that the analysis of stability and chaos (i.e. the “dynamics”) of real polynomial maps with real fixed points is reduced to the study of the canonical polynomial maps defined above, since we can always take any polynomial in to its canonical form by means of , determine the stability properties and then go back to the original polynomial. A commutative diagram of the conjugacy is
| (5.11) |
5.1 Stability and chaos in the canonical map of degree n
The derivative of , recalling to simplify notation, is
| (5.12) |
Evaluating (5.12) in the fixed point we get the eigenvalue function for each
| (5.13) | ||||
Then, the asymptotic stability condition implies that
| (5.14) |
From (5.14) we can recover all the stability conditions for the fixed points of the Canonical Quadratic and Cubic Maps.
which leads us to the following
Definition 5.5 (Product Distance Function).
Let be the Canonical Polynomial Map of -th degree and its fixed points, all of which depend upon the parameter . Let be a real fixed point among the latter. Then
| (5.15) |
is called the Product Distance Function (PDF) of .
The definition is motivated by the fact that is a product of the distances to of each of the other fixed points and that this quantity is fundamental in determining the stability of the fixed points. These distances are positive when and is negative when . We have stressed the dependence on the parameter in the definition of so that its character as a function is clear, stemming from the corresponding dependence on of the fixed points. In this way, the stability condition for the -th fixed point is reduced to
| (5.16) |
Since must be negative in order for to be stable as a sufficient condition, and an odd number of factors must be negative for the product in to be negative, it follows that, if is even or , an odd number of negative factors is a sufficient condition for the hyperbolic fixed point to be stable; i.e. if , an odd number of fixed points must lie below and, consequently, an even or zero (respectively, odd) number of fixed points must lie above if is even (respectively, odd). By similar arguments, we can prove
Proposition 5.2 (Necessary conditions for the stability of ).
Let , be defined as above and be a hyperbolic real fixed point of . The following are necessary conditions for to be an asymptotically stable fixed point:
-
•
if or is even, an odd number of fixed points must have values lower than ; or
-
•
if is odd and , zero or an even number of fixed points must have values lower than .
We must remark that the above conditions are not sufficient for a fixed point to be an attractor. The sufficient condition, however, is stated as
Theorem 5.3 (Sufficient condition for the stability of ).
Let , be defined as above and be a hyperbolic real fixed point of . Then, a necessary and sufficient condition for to be an attractor is that
| (5.17) |
Below the value of -2 there are other “stability bands” that lead to further period doubling bifurcations of the fixed points as they are crossed, but they must be calculated numerically and, as we have seen, depend on the degree of the polynomial.
5.2 Example
We will now give an example in which the above theory is applied to different degree polynomial maps.
Example 5.1 (Quartic Maps).
Using definition 5.2 for , we have that
| (5.18) |
Suppose has at least one real fixed point. Without loss of generality, suppose this fixed point is . Then
Making the substitution we can verify that we get
where
The stability a real fixed point is given by the product distance function
whose value must remain between minus two and zero in order for to be asymptotically stable; that is, if all fixed points are real,
| (5.19) | |||||
for 0, and to be asymptotically stable fixed points, respectively.
For example, let
| (5.20) | ||||
The plots of these fixed points with their corresponding parametric dependence on are shown in the lower panel of figure 5.1.
In light of proposition 5.2 we expect only and to be able to be asymptotically stable fixed points in any given range of . As the middle panel of figure 5.1 shows, precisely and are the fixed points whose product distances cross the stability band in the range of being plotted. As we recall, the product distance functions are the “stability conditions” of the fixed points. As long as the product distances remain within the stability interval, the fixed points are attractors, as we can verify in the upper panel of figure 5.1; also in this last panel, we can see the two attracting fixed points at the beginning of the plotted range; then, first loses its stability and gives rise to the period doubling bifurcations cascade which leads to chaotic behavior; later, zero also loses its stability and also gives rise to period doubling and chaos.
Chapter 6 An example application of the canonical map theory
We shall now consider an example application of the theory developed in chapters 3 through 5 to a couple of well known examples of discrete dynamical systems. This applications demonstrate a new method of solution for analyzing the dynamics of real polynomial maps that simplifies upon the procedures found in current literature.
6.1 The logistic map
The logistic map is the most immediate and obvious example application. We had already encountered briefly the logistic map, , in subsection 1.8 of chapter 1 and again in example 2.1 of chapter 2, where we saw that the map undergoes a series of period doubling bifurcations starting at the value of , ultimately achieving a chaotic nature at [Elaydi(2008), p. 47]. We will now proceed to utilize the theory developed in the previous chapters, particularly in chapter 3 and generalized in chapter 5, to this particular map.
We already know (or can easily calculate) that the fixed points of the logistic map are and [Elaydi(2008), p. 43]. With some minor algebra, we can then check that in the corresponding Linear Factors Form of the logistic map is then
| (6.1) |
| (6.2) | ||||
| (6.3) | ||||
| (6.4) | ||||
| (6.5) |
The corresponding non-zero fixed point of the canonical logistic map, defined in (3.12), is then simply
| (6.6) |
from which we can state that the canonical logistic map takes the explicit form
| (6.7) |
In order to determine the stability properties of these fixed points, both zero and nonzero, in the canonical logistic map, it is then sufficient, as we have proved in chapters 3 and 5, to observe the behavior of the Product Distance Functions (PDFs) of these fixed points, viz.
| (6.8) | ||||
| (6.9) |
By determining when these PDFs cross the stability bands whose boundaries are shown in table 3.1 we can readily determine when these fixed points are stable or unstable, when they bifurcate and when they reach any attracting periodic orbit for any , up to crossing the band. This whole process is depicted in figure 6.1. In particular, we can see from table 3.2, and again at figure 6.1, that when , the zero fixed point is attracting since its PDF lies within the first stability band, then exchanges stability at , when this last FP becomes stable and proceeds to a period doubling cascade upon its PDF, , crossing the bands defined by the bifurcation values , , etc. until reaching . This last value, by the way, agrees quite well and improves upon the approximation reported in [Elaydi(2008)] of 3.570 for the logistic map, since with the calculations of the present work . Finally, since these maps are topologically conjugate and it is known that the logistic map is chaotic starting with [Elaydi(2008)], we conclude that the CQM must be so starting from , which we may define as , our final “bifurcation” value.
6.2 Harvesting strategies
The connection of the logistic map with population models is old and well known. In the book by [Sandefur(1990)], there are a few examples of uses of second degree polynomials as recurrence functions which are used to model “harvesting” –or hunting– strategies of a population of animals. The main idea is that the population of animals grows whenever there is food and resources in the environment which, by account of its finite resources has a certain “carrying capacity”; this leads to a maximum population this environment can hold, the population growing according to the logistic model. The population may then be “harvested” –or hunted– at a certain rate yet to be specified and, depending on this rate, it is not hard to imagine that the final fate of the population of animals may be (i) extinction, if the rate is too high; (ii) steady population below the carrying capacity, if the rate is “just right”; or (iii) steady population at its maximum value dictated by the carrying capacity of the environment. A final fourth possibility –perhaps more rare– is to (iv) bring more individuals of the species from outside the system under analysis and introduce them to it, therefore making it possible for the population to surpass in number the carrying capacity of the system, but only to return naturally to the maximum value after a finite number “time-steps”. To examine this in detail consider the system defined by the recurrence relation
| (6.10) |
Here the population growth in any given period is proportional to both the initial population and the difference between this and the maximum population, normalized to the value of 1. The proportion constant is the growth rate . After some simplification, we can rewrite system as
| (6.11) |
We have yet to add the “harvesting term(s)”, which we might do so in several ways.
If we consider a fixed rate, say each period we harvest a proportion of the population (in terms of the maximum) then
| (6.12) |
is the corresponding map of this discrete dynamical system. From here we see that the Fixed Points Polynomial (FPP) of is111Remember for the -th degree polynomial.
whose FPs we can determine to be
from which we immediately calculate the Linear Factors Form (LFF) to be
| (6.13) |
From this form it is also straightforward to determine that in the corresponding canonical form we have
Both systems are related then by the linear transformations and . With this choice, the zero fixed point in the canonical map corresponds to and the nonzero fixed point to . To analyze the stability of this system, consider to be fixed and given and to be the parameter of this family of systems. The one immediate conclusion is that, for , we necessarily have that . Now, real implies . Over the value becomes complex which does not give any fixed points (but would mean “over-harvesting”). Since “harvesting” cannot be negative, it is clear corresponds to the maximum value of and that when . Remember now that to analyze the stability of we consider its PDF, . Analogously, . Since the maximum value of is and, being this the “unconstrained growth rate”, , we have that , therefore putting the in the stability range between and determined by the first stability band. This range is never left in any situation with physical meaning so we conclude that the nonzero fixed point, i.e. in the original system, is always the only asymptotically stable fixed point and the zero fixed point, i.e. in the original system, is always unstable. The only case to analyze with care is since then the two fixed points collide, but proposition 3.9 (on page 3.9) guarantees that in this case is semistable from the right.
In conclusion:
-
1.
When the population faces extinction asymptotically. Over this value extinction is achieved in a finite number of steps, there not being any more fixed points.
-
2.
When , and the population tends to .
-
3.
When , i.e. no harvesting, still tends to from the right and, correspondingly, the population tends asymptotically to .
6.3 Remarks on the application of the method
As somewhat instructive as the above examples may have been, the true simplifying power of the Canonical Map Theory is not yet fully met. Therefore we encourage the reader to apply this method to test the dynamics of any discrete dynamical systems family with a polynomial map depending on a single parameter, and evaluate the usefulness of the method by yourself. Families of maps with particular parametric dependencies and of cubic order will give the most useful applications for the time being, since the theory of higher degree polynomials is so far restricted to analyzing the stability of fixed points only.
Chapter 7 Conclusions and future work
Finally, lets sum up the main findings of this work and give some directions on the future work that arises from it. Basically, we can summarize the findings of this work as having successfully given conditions for the stability of the fixed points of any real polynomial map with real fixed points and that depends on a single parameter. In order to do this we have defined “canonical polynomial maps” which are topologically conjugate to any polynomial map of the same degree with real fixed points. Then, the stability of the fixed points of the canonical polynomial maps has been found to depend solely on a special function called “product distance” of a given fixed point. The values of this product distance determine the stability of the fixed point and when it bifurcates to give rise to attracting periodic orbits of period for all , and even ultimately when chaos arises through the period doubling cascade, as it passes through different “stability bands”, although the exact values and widths of these stability bands are yet to be calculated for regions of type greater than one for higher order polynomials. The latter must be done numerically.
As specific goals achieved we can remark that in this work we have
-
•
Generalized the existing results on parametric dependency of the coefficients of real quadratic maps for a larger set of functions.
-
•
Studied in detail the dynamics of general real cubic maps when the parametric dependency of its coefficients is given by continuous functions.
-
•
We have obtained transformations that take general quadratic and cubic maps into one whose dynamics is known and easy to study, and that are also topologically conjugate maps, so that the properties of stability and chaos are preserved.
-
•
Explicitly given the set of fixed points for the the proposed ’canonical’ quadratic and cubic maps.
-
•
Created a method that allows to understand the dynamics of cubic maps, that also includes quadratic maps.
-
•
Generalized the latter method to -th degree real polynomial maps.
-
•
Proposed a method to generate a large class of regular-reversal cubic maps.
-
•
Reproduced the know stability characteristics of the logistic map through the proposed method in a simpler way.
-
•
Proposed a methodology that allows to be used to create discrete dynamical systems with some prescribed bifurcation diagram.
As future work we can point out the following:
-
•
The numerical values of the limits of the stability bands for still higher degree canonical polynomial maps must be calculated with numerical methods.
-
•
Ultimately it is desired to obtain extensive tables of the bifurcation values for higher order polynomials.
-
•
Algebraic multiplicity of the roots of the Fixed Points Polynomial, when nonhyperbolic fixed points take place, must be analyzed for the general n-th degree Canonical Polynomial Map.
-
•
The case of complex fixed points must be included in detail for completeness.
-
•
Once the preceding is accomplished, the possibility of extending the theory to apply it to the Taylor Polynomial of any real analytical iteration function can be analyzed.
-
•
The power and simplicity of the proposed methodology will best be appreciated with 3rd or higher degree polynomials and when the implications on the Taylor polynomial of any non-linear map are understood.
Bibliography
- [Banks et al.(1992)Banks, Brooks, Cairns, Davis, and Stacey] J. Banks, J. Brooks, G. Cairns, G. Davis, and P. Stacey. On Devaney’s definition of Chaos. Am. Math. Month., 99:332–334, 1992.
- [Devaney(1989)] Robert Devaney. Introduction to Chaotic Dynamical Systems. Addison-Wesley Publishing Company, Inc., 1989.
- [Elaydi(2000)] Saber N. Elaydi. Discrete Chaos. Chapman & Hall/CRC, 2000.
- [Elaydi(2008)] Saber N. Elaydi. Discrete Chaos: with Applications in Science and Engineering. Chapman & Hall/CRC, Taylor & Francis Group, 2nd. edition, 2008.
- [Feigenbaum(1978)] Michael J. Feigenbaum. Quantitative Universality for a Class of Nonlinear Transformations. Journal of Statistical Physics, 1978.
- [Holmgren(1994)] Richard A. Holmgren. A First Course in Discrete Dynamical Systems. Springer-Verlag New York, Inc., 1994.
- [Marsden(1974)] Jerrold E. Marsden. Elementary Classical Analysis. W. H. Freeman and Company, San Francisco, CA, USA, 1974.
- [Protter(1998)] Murray H. Protter. Basic Elements of Real Analysis. Springer-Verlag New York, Inc., 175 Fifth Avenue, New York, NY 10010, USA, 1998. ISBN 0-387-98479-8.
- [Sandefur(1990)] James T. Sandefur. Discrete Dynamical Systems: Theory and Applications. Oxford University Press, 1990.
- [Singer(1978)] David Singer. Stable orbits and bifurcation of maps of the interval. SIAM J. Appl. Math., 2(35):260–267, 1978.
- [Solís and Jódar(2004)] Francisco J. Solís and Lucas Jódar. Quadratic Regular Reversal Maps. Discrete Dynamics in Nature and Society, 2004(2):315–323, 2004.
- [Spiegel et al.(2009)Spiegel, Lipschutz, and Liu] Murray R. Spiegel, Seymour Lipschutz, and John Liu. Mathematical Handbook of Formulas and Tables. McGraw-Hill, third edition, 2009. doi: 10.1036/0071548556.
- [Vellekoop and Berglund(1994)] Michel Vellekoop and Raoul Berglund. On Intervals, Transitivity = Chaos. Am. Math. Month., 101(4):353–355, 1994. URL http://links.jstor.org/sici?sici=0002-9890%28199404%29101%3A4%3C353%3AOIT%3DC%3E2.0.CO%3B2-B.
Appendix A Computational and numerical tools
For revision of the algebraic calculations, scripts in the wxMaxima scripting language were used, running under Ubuntu GNU/Linux 11.10. Both wxMaxima and Ubuntu are free software111“Free” as in “free speech”.. wxMaxima is a full-featured Computer Algebra System (CAS) in which some of the plots for this work were also made. Also, for most of the plots and bifurcation diagrams, a Python script shown below was used to calculate and plot most of the bifurcation diagrams presented in this work. The Python script was also ran under Ubuntu GNU/Linux and it uses numpy and matplotlib as libraries for numerical calculations and plotting, respectively.
import numpy as np
import matplotlib.pyplot as plt
import time
###################### DEFINITION OF FUNCTIONS #################################
def gn(n, xi, s=1, x0=0.1, niter=512, nplot=256): #iterates the canonical cubic map
#niter: total number of iterations
#nplot: number of points to return (for graphing)
#x0: initial condition
#xi: array containing the 3 fixed points
#n: order of the canonical map and number of f.p.’s
#s: sgn(M) from the original polynomial map
x = np.zeros(nplot) #here we store the orbit to be returned for graphing
xtemp = x0 #initialize orbit with initial condition
ssn=np.power(-1,n-1)*np.power(s,n) #the sign parameter
#part of the orbit that will not be graphed:
for i in range(niter-nplot):
prod=1
for j in range(n-1):
prod*=xtemp-xi[j+1]
xtemp *= (1 + ssn*prod)
x[0]=xtemp
#part of the orbit that will be graphed:
for i in range(nplot-1):
prod=1
for j in range(n-1):
prod*=x[i]-xi[j+1]
x[i+1] = x[i]*(1 + ssn*prod)
return x
def pdfn(n,xi): #Product Distance Function Values Array
#xi: vector containing the n fixed points
#k: number of fixed point whose stability we want to calculate
#n: total number of fixed points
pdf=np.zeros(n) #array to store product distance values for each f.p.
for k in range(n):
prod=1 #initalize the product
for j in range(n): #perform product
if(j != k):
prod *= xi[j] - xi[k]
pdf[k]=prod
return pdf
def xl(i, p, lamb): #define the fixed points as a function of the parameter
#i: number of fixed point to return value
#p: NP 2-d array containing the extra parameters for the fuctional forms
# of the fixed points. 1st index is the number of fixed point, 2nd is
# the number of relevant parameter for that fixed point
if(i>0):
#return lamb*(3-lamb) #quadratic FP for QCM
#return p[i][1]*np.sqrt(lamb) + p[i][2] #square root fixed points
return p[i][0]*lamb + p[i][1] #linear fixed points
#return (-1)**p[i,0]*p[i,1]+(-1)**(p[i,0]+1)*np.exp(-p[i,2]*lamb)
else:
return 0 #the zeroth fixed point is always zero in value
####################### DEFINITION OF PARAMETERS ###############################
print time.clock(), "Defining parameters..."
n=4 #number of fixed points and order of the canonical map
npl=256 #number of points to graph for each value of lambda
nit=512 + npl #total number of iterations
lmin=0 #starting value of lambda for the graph
lmax=1.2 #ending value of lambda for the graph
p=np.zeros([n,3]) #parameters for the definition of the fixed points
p[1,0]=1
p[1,1]=0
#p[1,2]=-1
p[2,0]=-1
p[2,1]=0
#p[2,2]=1
p[3,0]=2
p[3,1]=0
#p=np.array([[0,0,0],[1,0,0],[-1,0,0]]) #for the linear FPs
#p=np.array([[0,0,0],[1,0,0],[-1,0,0]]) #for the linear FPs
#p=np.array([[0,0,0],[0,2,0],[6,1,0]]) #one linear FP another const. FP
lval=np.linspace(lmin,lmax,600) #range of the parameter for the bif. diagram
midl=(lmin+lmax)/2 #middle point of the lambda interval (bif.val. estimation)
uncert=abs(lmax-lmin)/2 #uncertainty for the bif. value estimation
xi=np.zeros(n) #the zero fixed point is always zero in value
#x0 = 0.05 #initial condition
pdf=np.zeros([n,lval.size])
fps=np.zeros([n,lval.size])
lines=[’-’,’--’,’-.’]*3
linw=[2,2,2,3,3,3,4,4,4,5,5,5]
############# PRINT ESTIMATIONS TO SCREEN ######################################
print time.clock(), "Print parameter estimates..."
#we print the values found for the bifurcation
#along with its corresponding propagated uncertainties
print "l_i \in (",lmin,", ",lmax,")"
print "l_i =",midl,"+-",uncert
#print "b_i \in (", -p[1][1]-np.exp(-lmin),",",-p[1][1]-np.exp(-lmax),")"
#print "b_i =",-p[1][1]-np.exp(-midl),"+-",midl*np.exp(-midl)*uncert
############# BEGIN PLOT #######################################################
print time.clock(), "Beginning Fixed Points plot..."
##draw the graph of the n fixed points varying with the parameter
plt.figure(1)
plt.subplot(313)
plt.grid(True)
plt.xlabel(’Parameter’)
plt.ylabel(’FPs’)
plt.xlim(lmin,lmax)
#plt.title(’Fixed Points’)
#plt.ylim(-2,2)
for k in range(n):
fps[k]=[xl(k,p,lam) for lam in lval]
plt.plot(lval,fps[k,:],’k’,label=’FP’+str(k), ls=lines[k],linewidth=linw[k])
plt.legend()
print time.clock(), "Beginning PDFs plot..."
#draw the product distance functions of each of the n FP’s
plt.subplot(312)
plt.grid(True)
#plt.xlabel(’Parameter’)
#plt.ylabel(’PDFs’)6
#plt.ylabel(’Stability Conditions’)
plt.ylabel(’PDFs’)
#plt.title(’PDFs of the Fixed Points’)
plt.xlim(lmin,lmax)
plt.ylim(-4,1)
#print bifurcation values boundaries
plt.plot(lval,np.zeros_like(lval),’k--’)
plt.plot(lval,-2*np.ones_like(lval),’k--’)
if(n==2):
plt.plot(lval,-2.449*np.ones_like(lval),’k--’)
plt.plot(lval,-2.544*np.ones_like(lval),’k--’)
plt.plot(lval,-2.5642*np.ones_like(lval),’k--’)
#plt.plot(lval,-2.56871*np.ones_like(lval),’k--’)
#plt.plot(lval,-2.56966*np.ones_like(lval),’k--’)
#plt.plot(lval,-2.569881*np.ones_like(lval),’k--’)
plt.plot(lval,-2.57*np.ones_like(lval),’k--’)
if(n==3):
plt.plot(lval,-3*np.ones_like(lval),’k--’)
plt.plot(lval,-3.236*np.ones_like(lval),’k--’)
plt.plot(lval,-3.288*np.ones_like(lval),’k--’)
plt.plot(lval,-3.29925*np.ones_like(lval),’k--’)
for j in range(lval.size):
pdf[:,j]=pdfn(n,fps[:,j])
for k in range(n):
plt.plot(lval,pdf[k,:],’k’,label=’PDF’+str(k),ls=lines[k],linewidth=linw[k])
#plt.plot(lval,pdf[k,:],’k’,label=’SC’+str(k),ls=lines[k],linewidth=linw[k])
plt.legend()
print time.clock(), "Begin bifurcation diagram plotting..."
#draw the bifurcation diagram
plt.subplot(311)
plt.grid(True)
#plt.xlabel(’Parameter’)
plt.ylabel(’Asymptotic Value’)
plt.title(’Bifurcation Diagram’)
plt.xlim(lmin,lmax)
#plt.ylim(-2, 2)
m=0
for lamb in lval: #draws the bifurcation diagram
# for j in range(n-1): #calculate the fixed points
# xi[j+1]= xl(j+1,p,lamb)
#
# aux=pdfn(n,xi) #calculate PDF of the FP’s
# for k in range(n):
# pdf[k][m]=aux[k] #store the values for graphing
#plt.subplot(311)
x = gn(n, fps[:,m], x0=0.9, niter=nit, nplot=npl)
plt.plot(lamb*np.ones(npl),x,’k,’)
x = gn(n, fps[:,m], x0=-0.01,niter=nit, nplot=npl)
plt.plot(lamb*np.ones(npl),x,’k,’)
m+=1
print time.clock(), "Done."
plt.show()